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Peptide-MHC (pMHC) ligand engagement by T-cell receptors (TCRs) elicits a variety of cel- 
lular responses, some of which require substantially moreTCR-mediated stimulation than 
others. This threshold hierarchy could reside at the receptor level, where different response 
pathways branch off at different stages of theTCR/CD3 triggering cascade, or at the cellu- 
lar level, where the cumulative TCR signal registered by the T-cell is compared to different 
threshold values. Alternatively, dual-level thresholds could exist. In this study, we show 
that the cellular hypothesis provides the most parsimonious explanation consistent with 
data obtained from an in-depth analysis of distinct functional responses elicited in a clonal 
T-cell system by a spectrum of biophysically defined altered peptide ligands across a range 
of concentrations. Further, we derive a mathematical model that describes how ligand 
density, affinity, and off-rate all affect signaling in distinct ways. However, under the kinetic 
regime prevailing in the experiments reported here, the TCR/pMHC class I (pMHCI) dis- 
sociation rate was found to be the main governing factor. The CD8 coreceptor modulated 
the TCR/pMHCI interaction and altered peptide ligand potency. Collectively, these find- 
ings elucidate the relationship betweenTCR/pMHCI kinetics and cellular function, thereby 
providing an integrated mechanistic understanding of T-cell response profiles. 
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INTRODUCTION 

T lymphocyte antigen receptors mediate adaptive immune 
responses via interactions with disease-associated peptide ligands 
presented on the surface of target cells by major histocompatibility 
complex (MHC) molecules. In the case of CD8 + T-lymphocytes 
(CTLs), which constitute the classical T-cell effector subset, the 
clonotypically expressed T-cell receptors (TCRs) engage spe- 
cific peptide-MHC class I (pMHCI) molecules to elicit several 
functions that are instrumental in eliminating the pathogenic 
threat (1). 

The six hyper-variable complementarity-determining regions 
(CDRs) of the TCR govern molecular interactions with pMHCI 
(2). These CDRs confer the specificity of molecular recognition, 
allowing CTLs to attack diseased cells without causing undue harm 
to healthy cells (1). Nonetheless, a certain degree of degeneracy 
within the antigen recognition system is unavoidable (3). Indeed, 
ample experimental evidence supports the notion that a single 
TCR clonotype can interact productively with numerous peptide 
ligands, which typically vary in their ability to elicit different types 
of cellular response (4-14). Furthermore, some cellular responses 
are more readily evoked than others (15-18). It seems reasonable 
to explain this phenomenon by postulating a hierarchy of thresh- 
olds. Such a hierarchy may reside either at the level of individual 
TCR/CD3 complexes or at the cellular level. In the former sce- 
nario, the threshold parameters relate to distinct components of 



the TCR triggering process, with different responses being elicited 
at different stages in the development of a mature signalosome. 
Essentially, this is the kinetic discrimination model proposed by 
Rabinowitz et al. (19), whereby early responses require a shorter 
TCR/pMHCI dwell-time than late responses. This model is logi- 
cally distinct from the kinetic proofreading model that accounts 
for the existence of dwell-time thresholds per se (20). 

In the cellular-level scenario, TCR triggering delivers a stereo- 
typical signal that elicits distinct responses across a single quantita- 
tive gradient (e.g., concentration of the relevant signaling factor). 
The contrast between these two scenarios is illustrated diagram- 
matically in Figure 1. A third possibility is that the hierarchy 
comprises a combination of both receptor-level and cellular- 
level modes. It has hitherto been unclear which of these three 
alternatives prevails in T-cell signaling. 

In this study, we used a mathematical model to investigate the 
extent to which each of the three possibilities (receptor-level hier- 
archy, cellular-level hierarchy, or a combination of both) agree 
with experimental data. Our model accounts for the functional 
sensitivity of TCR-mediated responses to a given pMHCI ligand 
on the basis of TCR/pMHCI interaction kinetics. The concept 
that TCR/pMHCI kinetics governs functional sensitivity was pio- 
neered by Lanzavecchia et al. (21) and subsequently formulated 
as a mathematical model (22, 23). Bridgeman et al. (24) recently 
summarized the available evidence across reported systems. 
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FIGURE 1 | Two hypotheses on the level of the response hierarchy. (A) 

Signaling cascades are initiated at different points along theTCR/CD3 reaction 
pathway, which includes phosphorylation events and the docking of kinases 
and linker proteins. Thus, responses R u R 2 , and R 3 would correspond to 
distinct values of theTCR triggering threshold, i.e., F™ < Tj 2 ' < Tj 31 . (B) 
Cellular integration of signals from the triggeredTCRs is shown as a 
summation (£) box. This signal passes through a non-linear threshold that 



determines whether the cell will respond. Here, R u R 2 , and ff 3 represent 
various cellular responses, such as the expression of different cytokine species. 
Each receives the same integrated signal Was an input, but different values of 
Ware required to initiate a response. This is depicted schematically by sigmoid 
curves whose midpoint lies to the left, in the middle, or to the right. These 
differences correspond to distinct values of the cellular activation threshold, 
i.e., Will < Will < H/, 13 , 1 . These two hypotheses are not mutually exclusive. 



Analysis of TCR/pMHCI interaction kinetics demonstrates that 
functional sensitivity is not dependent on a single biophysical 
parameter, but rather on the interplay between association rate, 
dissociation rate, and ligand densities (23, 25, 26). Under the 
experimental conditions prevailing in the present study, however, 
the dissociation rate (the reciprocal of the mean dwell-time of 
the TCR/pMHCI interaction) emerged as the dominant biophys- 
ical parameter. Moreover, we estimated the extent to which the 
CD8 coreceptor modulates this parameter using data from par- 
allel experiments conducted in the absence of an extracellular 
MHCI/CD8 interaction. 

Five distinct cellular responses were investigated: (i) mobi- 
lization of the degranulation marker CD107a; (ii) secretion of 
macrophage inflammatory protein 1-/3 (MIP-1/3); (iii) secretion 
of tumor necrosis factor-a (TNF-a); (iv) secretion of interleukin- 
2 (IL-2); and (v) secretion of interferon-y (IFN-y). Measured 
simultaneously and independently by flow cytometry, the magni- 
tude of each response was quantified as a function of fluorescence 
and plotted against peptide concentration. The resulting dose- 
response curves indicated that the five cellular responses were 
elicited by pMHCI stimulation according to a pronounced hierar- 
chy. Analysis of these curves by means of the mathematical model, 
in conjunction with biophysical data, indicated that the cellular- 
level response threshold hierarchy hypothesis provides the most 
parsimonious explanation. 

MATERIALS AND METHODS 
EXPERIMENTAL PROCEDURES 
Cells and peptide ligands 

The CTL clone ILA1 recognizes residues 540-548 (ILAKFLHWL 
in single-letter amino acid code, abbreviated as ILA) of the human 
telomerase reverse transcriptase protein presented in the context of 
the human MHCI allotype HLA-A*0201 (HLA-A2). Cell culture 



was performed as described previously (27). The altered peptide 
ligands used in this study are referred to here as 3G, 5Y, 8T, 8E, and 
3G8T; these ligands were largely characterized previously (27, 28), 
and display equivalent binding to HLA-A2 (29). 

Bioassay 

Stable C1R transfectants expressing wild-type HLA-A2 (C1R- 
A2) or CD8-null HLA-A2 (C1R-A2null), the a3 domain of 
which contains the double mutation D227K/T228A that abro- 
gates CD8 binding (30), were pulsed with peptide as indi- 
cated for 1 h at 37°C. For each condition, assays were set up 
simultaneously in 96-well plates using 4.5 x 10 5 C1R cells per 
well, thereby ensuring that all cellular parameters were con- 
sistent across ligand stimulations. Subsequently, brefeldin A 
(10/xg/ml; BD Biosciences) and monensin (0.7 /xl/ml; BD Bio- 
sciences) were added together with a pre-titred concentration of 
the directly conjugated monoclonal antibody (mAb) aCD107a- 
FITC [BD Pharmingen; Ref. (31)]. Serum-starved CTLs, incu- 
bated in Roswell Park Memorial Institute (RPMI) medium (Life 
Technologies) containing 2% fetal calf serum (FCS) for 16 h 
prior to assay, were subsequently added at 9 x 10 5 cells per 
well. After incubation for 6h, which generally allows sufficient 
time for even the most sluggish response to appear, cells were 
washed in phosphate-buffered saline (PBS) containing 1% FCS 
and incubated for lOmin at room temperature with Aqua fix- 
able live/dead cell stain (Life Technologies) to enable the exclu- 
sion of dead cells from the analysis. Pre-titred concentrations of 
aCD3-H7allophycocyanin (BD Biosciences), aCD8-QD705 (Life 
Technologies), and aCD19-V500 (BD Horizon) were then added 
for 20min at 4°C. After two further washes in PBS/1% FCS, 
cells were fixed/permeabilized using a Cytofix/Cytoperm™ kit 
(BD Biosciences) according to the manufacturer's instructions, 
then stained intracellularly with pre-titred concentrations of 
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aMIP-l^-PE (BD Pharmingen), aIFNy-V450 (BD Horizon), 
aIL-2-allophycocyanin (BD Pharmingen), and aTNFa-PECy7 
(BD Pharmingen) for 20 min at 4°C. 

Flow cytometry 

Stained cell samples were acquired and recorded using a cus- 
tomized FACS Aria II flow cytometer (BD Biosciences) equipped 
for the simultaneous detection of 18 fluorescent parameters. Poly- 
chromatic analyses were conducted using Flowjo software version 
9.5.2 (TreeStar Inc.). The following gating tree was applied: (i) 
single cells were identified based on their light scatter proper- 
ties; (ii) Boolean gating was carried out to exclude artifacts and 
fiuorochrome aggregates; (iii) viable CD3+CD19~ events were 
selected; (iv) outliers were eliminated in a side-scatter versus CD3 
display; and (v) gates were set on cells positive for individual 
functional read-outs. The frequencies and median fluorescence 
intensity (MFI) values of each functional response were exported 
for data analysis to Excel and onward for simultaneous curve fitting 
in Mathematica. 

Surface plasmon resonance 

Soluble TCR, derived from the ILA1 CTL clone, was manufac- 
tured as described previously (32, 33). Binding analysis by surface 
plasmon resonance (SPR) was performed using a BIAcore T100™ 
equipped with a CM5 sensor chip (34). Biotinylated pMHCI (200- 
400 response units) was immobilized to streptavidin, which was 
chemically linked to the chip surface. The pMHCI was injected at 
a slow flow rate ( 1 0 /A/mm) to ensure uniform distribution on the 
chip surface. Combined with the small amount of pMHCI bound 
to the chip surface, this reduced the likelihood of dissociation rate 
limiting mass transfer effects. The ILA1 TCR was purified and con- 
centrated to ~100/iM on the day of SPR analysis to reduce the 
likelihood of TCR aggregation affecting the results. At least 5 serial 
dilutions were prepared in duplicate and injected over the relevant 
sensor chips at a flow rate of 45 fil/min. All experiments were con- 
ducted at 25°C. Results were analyzed using BIAevaluation, Excel, 
and Origin. 

THEORY AND DATA ANALYSIS 
Assumptions 

The assumptions of the TCR triggering model are as follows: (i) 
TCRs on the T-cell surface become "triggered" (i.e., are induced to 
become signalosomes) during an interaction with a pMHCI lig- 
and; and (ii) the T-cell accumulates the signals emanating from 
triggered TCRs over space (celhcell interface) and time (dura- 
tion of cellxell interaction), proceeding with a response when this 
signal exceeds a cellular activation threshold. The mathematical 
formulation of these assumptions is based on the kinetics of the 
TCR/pMHCI interaction. 

The TCR triggering rate equation 

The model describes the kinetics of interactions between TCR 
and pMHCI molecules in the interface area between a T-cell and 
an antigen-presenting cell (APC). This area is occupied by TCRs 
on the T-cell side of the interaction and by pMHCI complexes 
on the APC side. Let Rj be the total number of TCRs and Zj 



be the total number of pMHCI molecules. Both are subject to a 
conservation law: 



R T = C + R 
Z T = C + Z 



(1) 
(2) 



where C is the number of TCR/pMHCI complexes, _R is the 
number of TCRs not engaged in a complex (i.e., "free" TCR mol- 
ecules) and Z is the number of free pMHCI molecules. Kinetic 
equilibrium is expressed by the law of mass action equation: 



k on RZ = k oS C. 



(3) 



This equation, together with equations (1) and (2), leads to a 
quadratic in C with solution 



C: 



R T + Zj + koff/kon 



4R T Z T 



R T + Z T + k oS /k 0 



(4) 

(the other root is irrelevant since it exceeds both Rj and Zj). The 
rate at which TCRs are being triggered can be expressed as the 
rate at which TCR/pMHCI complexes dissociate times the prob- 
ability that any given TCR/pMHCI association event results in a 
triggering of the associated TCR/CD3 complex: 



W : 



Koff 



CP, 



trig 



(5) 



where W denotes the TCR triggering rate and P tr ; g the triggering 
probability for an individual interaction event. Since triggering 
of a TCR/CD3 complex requires completion of a series of phos- 
phorylation and docking events, it is reasonable to assume that 
triggering can only happen if the peptide ligand remains engaged 
for at least a certain amount of time. This minimum duration is 
the TCR triggering threshold T R . If the lifetime of the TCR/pMHCI 
complex is exponentially distributed, the triggering probability is 
given by: 



P trig = expf-fcoffTR}. 



(6) 



Combined, equations (4-6) yield the general triggering rate 
equation: 

t . r ^off exp{-fcoff T R ] 

W = (R T + Z T + fc 0 ff/fcon) 



x 1-1 



R T + Z T + Ks/kon 



(7) 



Simple forms in special cases [equations (13-15)] are obtained 
by using the first-order approximation to the square root in 
equation (7). 

Connecting the theoretical model to experimental observations 

The available data comprise cellular response measurements 
(read-outs) of T-cell activation following exposure to APCs incu- 
bated with various peptide ligands across a range of concen- 
trations, as well as TCR/pMHCI association and dissociation 
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rates for each peptide ligand, measured via SPR. To relate these 
measurements to the TCR triggering theory, a number of auxiliary 
assumptions are required. In particular, if Y is the peptide incu- 
bation concentration, the following proportionality is assumed: 



Z T = aY. 



(8) 



If Ko is the dissociation constant determined via SPR, the 
following proportionality is assumed: 



fcoff 



(9) 



The proportionality constant k is required to convert rates 
as measured via SPR to the 2-dimensional environment of the 
cellxell interaction area. To connect the model to data, two com- 
pound parameters are introduced: f =a/K (dimensionless), and 
p = RjIk (M). With these parameters, the TCR triggering rate 
w = W/k (in M-s -1 ) assumes a simpler form: 

k oS exp{- fc off T R ] 
w = (p + f Y + K D ) 



x 1 



4p$Y 



p + ZY + K D 



(10) 



The read-out in the present study is median fluorescence inten- 
sity (MFI), which is a valid measure of activation in view of its 
good correlation with the fraction of responding cells. Accord- 
ingly, the read-out U is assumed to be proportional to the fraction 
of responding cells. This assumption is validated by response data 
at the level of individual cells (see Results below). The equation 
for the read-out is as follows: 



U = U min + (U a 



U min ) P(respond) 



(11) 



where U m i n and L7 max are nuisance parameters associated with 
the read-out procedure; these parameters take specific values for 
each type of cellular response and are assumed not to vary across 
peptide ligands. The probability that a T-cell responds is modeled 
as follows: 



P(respond) = P(w act < w) 



(12) 



where w is the scaled TCR triggering rate given by equation (10) 
and w act is the cellular activation threshold. This latter quantity is 
assumed to have a log-normal distribution over the population 
of T-cells in the experiment. Stochastic variation between the T- 
cells within the responding population is thus taken into account. 
The parameters are estimated by simultaneous least-squares fit- 
ting over the set of available ligands and read-outs. The median of 
the log-normal distribution provides the estimate of iv act . 

Implementation of response hierarchy hypotheses 

We examined five different cellular responses and compared three 
hypotheses: (i) for each type of cellular response, both Tr and 
w act have distinct values; (ii) there is a common value of Tr across 
responses, whereas w act has a different value for each response; and 



(iii) there is a common value of w act across responses, whereas 
Tr has a different value for each response. As the least-squares 
fit is carried out simultaneously, these hypotheses can be imple- 
mented by specifying either common or particular values for these 
parameters. 

RESULTS 

KINETIC PARAMETERS AND LIGAND NUMBERS INTERACT TO 
DETERMINE THE RATE OF TCR TRIGGERING 

The triggering equation (7) reduces to various simplified forms, 
depending on which receptors are kinetically limiting. These forms 
are of immunological interest since they illuminate the controversy 
over which kinetic parameter primarily governs T-cell activation. 
First, if the term k 0 f[/k on is much larger than both Rj and Zj, 
which applies when both cells have introduced comparatively 
low numbers of ligands into the interaction area, the following 
approximation is accurate: 

W = R T Z T k on exp{-k oS T R } (k oS /k on » max{_R r ,Z T }) , 

(13) 

that is, the triggering rate will be proportional to both TCR and 
pMHCI numbers: signaling is affinity-limited. 

On the other hand, if the cells express numbers of lig- 
ands in the interaction area that are much larger than the 2- 
dimensional dissociation constant k Q g/k on , there are two further 
good approximations: 

W = RtKs exp{-fcoff Tr} (koff/kon « Z T ) (14) 
W = Z T k off exp{-fc off T R } (fcoff/U « Rt); (15) 

in the first case, signaling is TCR-limited and in the second, it 
is MHC-limited. These equations imply that under the receptor- 
limited regimes, the TCR triggering rate displays an optimum 
relative to k 0 g at the point fc 0 g- = T^ 1 , whereas under the affinity- 
limited regime, W increases monotonically in k on and decreases 
monotonically in k 0 g . 

When both Rj and Zj are much larger than k 0 g/k on , the 
TCR triggering rate is proportional to min{J?j-, Zj} as shown in 
Figure 2A, whereas at lower receptor densities, W oc RjZj as per 
equation (13). The transition between the two regimes happens 
where the receptor densities traverse the 2-dimensional dissoci- 
ation constant k 0 qlk on . Thus, depending on receptor densities, 
the triggering rate may be affected solely by changes in Rj (TCR 
copy numbers), Zj- (pMHCI copy numbers), or both. The MHC- 
limited case corresponds to the serial triggering (serial engage- 
ment) hypothesis proposed by Valitutti et al. (35-37); however, 
the present theory is more general in the sense that the serial 
triggering mechanism arises as a special case. 

Figure 2B shows that the TCR triggering rate increases with 
both k on and mean dwell-time k~g when fc on is low (relative to 
receptor densities; scaling is explained in the legend to Figure 2). 
However, when fc on increases, the triggering rate exhibits a weak 
dependence on k on and a non-monotone dependence on k Q g, 
attaining a maximum at fc 0 g- = 1/Tr. 

The avidity effect is routinely exploited to determine the func- 
tional sensitivity of the TCR to a given ligand. This effect hinges 
on mutual compensation by pMHC ligand numbers and the TCR 
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FIGURE 2 |TCR triggering rate depends onTCR/pMHCI association and 
dissociation rates as well as the densities of TCR and pMHCI molecules 

(A) Dependence of the rate of TCR triggering on receptor densities. 
Quantities are dimensionless: TCR triggering rate is scaled as WT^k 0 „, TCR 
density is scaled as R T koJk M , pMHCI density is scaled as Z T k^/k m , and 
k M = 1/T B is assumed (this means that the ligand is optimal under 
MHC-limited conditions). (B) Dependence of the rate of TCR triggering on 
kinetic parameters. The dependence on the mean TCR/pMHCI interaction 
time (1/^ofi) is non-monotone, indicating that a maximally strong agonist 
must satisfy k M = 1/7",,. In contrast, the dependence on k m is monotone. 
Affinity is customarily expressed by the dissociation constant K D = k M /k„ . 
Accordingly, whether or not an improvement of affinity correlates with an 
enhanced TCR triggering rate depends on where the system is initially 
located on the surface of the graph, as well as on the relative contributions 
that changes in association and dissociation rates make to the overall change 
in affinity. Quantities are dimensionless: TCR triggering rate is scaled as 
WT„I(R T + Z T ), association rate is scaled as k m T R (R T + Z T ), and the mean 
lifetime of the interaction is scaled as T/(k M T R ). Here, R T =Z T . Similar 
qualitative behavior results when the two ligand densities are unequal. 



triggering rate per ligand; i.e., a T-cell can be activated by a strong 
agonist even at low copy numbers whereas a poor agonist may be 
effective but only at sufficiently high copy numbers. In equation 
(15), copy numbers are expressed by Zj and the intrinsic quality 
of the ligand is expressed by the quantity k 0 nexp{—k 0 f{Ts}. The 
latter attains a maximum at k Q g = 1/Tr. If the triggering rate W 
needs to exceed a certain threshold value to activate the T-cell, suf- 
ficient pMHCI copy numbers can compensate for a suboptimal 
intrinsic triggering rate. However, increasing Zj may transfer the 
system into the TCR-limited regime (Figure 2B) , under which con- 
ditions pMHCI density-based compensation cannot be achieved, 
cf. equation (14). Moreover, physiological bounds must exist on 
how many copies can be present in the interaction area and poor 



agonists will be unable to activate the T-cell unless unusually high 
numbers of the ligand are presented (e.g., the cell transcribes very 
large quantities of the protein). 

INTERCELLULAR VARIATION IN RESPONSE HIERARCHY THRESHOLDS 

An overview of the experimental data is provided in Figure 3. In 
keeping with the avidity effect, the fraction of cells that exhibited at 
least n distinct responses (where 1 < n < 5) increased with peptide 
concentration; individual ligands differed with respect to the min- 
imal concentration required to elicit n cellular responses. There 
was appreciable variability with regard to which functions were 
elicited, as depicted in Figure 3B, which displays the data obtained 
at peptide concentrations of 1CP 6 M. The responses generally 
conformed to the following series: 

MIP-l/S < CD107a < TNFa < IL-2 < IFN-y 

(with the lowest threshold on the left). However, a small propor- 
tion of cells exhibited a slightly different ordering of threshold 
values, as would be expected given the natural variability between 
cells with respect to the threshold value of any one given response. 

These results indicate that, for each response, there is varia- 
tion between cells with respect to the cellular threshold value. This 
supports the assumption made in the section Connecting the the- 
oretical model to experimental observations above, namely that the 
sigmoid dose-response shape of the population-level read-out can 
be accounted for by postulating a statistical distribution at the cell 
population level with respect to the cellular threshold for any given 
response. 

THE RESPONSE HIERARCHY RESIDES AT THE LEVEL OF CELLULAR 
SIGNAL INTEGRATION 

The parameters of the mathematical model were estimated by 
means of least-squares fitting of model predictions to the data; the 
biophysical data are summarized in Table 1 and Figure Al in the 
Appendix. Figure 4 shows the results of least-squares curve fit- 
ting for the dual-hierarchy hypothesis, the receptor-level response 
hierarchy hypothesis, and the cellular-level response hierarchy 
hypothesis. The goodness-of-fit should be assessed in the light 
of the parsimoniousness of the model. There is less than one para- 
meter per curve. This compares favorably to the standard practice 
of fitting a separate sigmoid model to each curve, which requires 
three or more parameters per curve. The fit is very good for such 
extreme parameter-count efficiency. 

There is no substantial difference in the quality of fit between 
the dual-threshold model and the cellular-threshold model, 
whereas the fit to the receptor- threshold model is consistently less 
good. Thus, the data are best explained by postulating a hierarchy 
only at the level of the cellular activation threshold. 

The MlP-l/i response was found to have the lowest acti- 
vation threshold. Relative to this value, the ratios of median 
threshold values were found to be as follows: CD107a/MIP- 
1/3 = 1.44; TNFa/MIP-1/3 = 11.7; IL-2/MIP-lyS = 225; and IFN- 
y /MIP- 1 p = 23 1 . These findings broadly agree with the consensus 
series derived from the data reported in Figure 3. The esti- 
mated median activation threshold values are markedly higher 
for IL-2 and IFN-y, which required a stimulus over two orders of 
magnitude larger compared to the more readily elicited responses. 
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FIGURE 3 | Functional variation at the cell population level in 
response to six different peptide ligands. The CTL clone ILA1 was 
stimulated for 6 h with peptide-pulsed C1 R-A2 targets as indicated. Five 
functional read-outs (CD107a, MIP-1,S,TNF-a, IL-2, and IFN-y) were 
measured via flow cytometry. The 8E peptide was included as an 



extremely weak ligand. (A) Overview of functional profiles. Pie chart 
segments represent the fraction of cells expressing the number of 
functions indicated in the key. (B) Detailed analysis of functional profiles at 
[peptide] = 10~ 6 M. Pie charts are extended with arcs defining expressed 
functions as indicated in the key. 



Table 1 | Biophysical parameters. 



Ligand 


On-rate (M- 1 s _1 ) 


Off-rate (s" 1 ) 


Kd (IV) 


3G8T 


19,500 


0.049 


2.5 x 10" 6 


3G 


16,000 


0.047 


2.9 x 10" 6 


8T 


4,000 


0.095 


2.4 x 10" 5 


I LA 


4,100 


0.13 


3.2 x 10" 5 


5Y 


1,300 


0.32 


2.5 x 10" 4 



The kinetics of the TCR/pMHCI interaction for the 5Y and 3G8T ligands are 
reported here for the first time; it should be noted that the 5Y data are subject to 
a comparatively greater potential measurement error due to their rapidity. 



MHCI/CD8 BINDING MODULATES FUNCTIONAL SENSITIVITY 

The experiments described above were repeated using CI R-A2null 
cells as APCs, which cannot productively engage the CD8 corecep- 
tor (38). It is well-attested that the MHCI/CD8 interaction mod- 
ulates the TCR/pMHCI association rate (k on ), the TCR/pMHCI 



dissociation rate (k Q fr ), and the receptor-level triggering threshold 
(Tr) (27, 38, 39). Multipliers can be used to model these effects. 
For example, k 0 f{ is replaced throughout by Koff^off where y Q ff rep- 
resents the effect of abrogating the MHCI/CD8 interaction and 
/off > 1 since CD8 stabilizes the TCR/pMHCI interaction (38- 
40). The parameter y Q ff was set to 2.16, the value reported by 
Wooldridge et al. (38). 

Since the prevailing kinetic regime in the experiments described 
here is MHC limitation, the parameter y on for the association rate 
could not be estimated and was set to 1. This leaves a single free 
parameter, y$_, for the data obtained with C1R-A2null APCs. Data 
fits are shown in Figure 5, with y& = 1.14 (the least-squares esti- 
mate) for the receptor threshold Tr. This value indicates that the 
MHCI/CD8 interaction decreases the TCR triggering threshold 
by 13%. 

DISCUSSION 

T-cells can exhibit a variety of cellular responses to TCR-mediated 
stimulation, some of which are elicited far more readily than 
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FIGURE 4 | Dose-response curves for five peptide ligands and five 
cellular responses. Solid black curves correspond to the dual-threshold 
model in which both the receptor-level and the cellular-level thresholds have 
different values for different cellular responses. Solid gray curves 
correspond to the receptor-level threshold model in which the receptor-level 



thresholds have different values for different cellular responses, but the 
cellular activation threshold is constant. Dashed curves correspond to the 
cellular-level threshold model in which the cellular activation thresholds have 
different values for different cellular responses, but the receptor-level 
threshold is constant. 




Peptide incubation concentration (Log 10 M) 



FIGURE 5 | Dose-response curves for five peptide ligands and five 
cellular responses in the absence of the MHCI/CD8 interaction All curves 
correspond to the cellular-level threshold model in which the cellular activation 
thresholds have different values for different cellular responses, but the 
receptor-level threshold is constant. Black curves are the fits to the data 



shown. Gray curves correspond to the situation in which the MHCI/CD8 
interaction is present. Parameter values are as in the previous figure. Two new 
parameters (ym and y K ) have been introduced and estimated to account for 
the absence of the CD8 effect; y M was set to 2.16 and y R was estimated as 
the sole remaining free parameter. 
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TCR/pMHCI dissociation rate (1/s) 

FIGURE 6 | Differential effects of CD8 modulation on functional 
sensitivity. TCR triggering rate, scaled with respect to the optimal ligand 
under MHC-limited conditions, as a function of k a „, comparing the wild-type 
MHCI/CD8 interaction ("A2") to the abrogated MHCI/CD8 interaction 
("A2null"). The positions for the five ligands are indicated. The inset shows 
two hypothetical ligands: "X," which is unaffected by the absence of the 
CD8 interaction, and "Y" which is optimal under the MHC-null condition. 



others. In particular, different responses require different levels of 
stimulation in terms of ligand densities present on the APC sur- 
face. Such threshold differences might be postulated at the receptor 
level (different response pathways branch off at different stages in 
the sequence of signaling events), and/or at the cellular level (the 
TCR signal received by the T-cell needs to exceed different thresh- 
old values for different responses). The present study suggests that 
the cellular-level hypothesis is the most parsimonious explanation 
consistent with the data. Furthermore, the role of the CD8 core- 
ceptor as a modulator of T-cell functional sensitivity is confirmed 
by the present findings, which indicate that the MHCI/CD8 inter- 
action decreases the TCR/pMHCI dissociation rate as well as the 
duration of TCR/pMHCI contact required to activate the T-cell. 

A key consideration for the interpretation of our results is the 
standardized nature of the system. Determinative cellular factors, 
such as TCR density, membrane constituency, and the expression 
of costimulatory/inhibitory molecules, were fixed across all condi- 
tions, which means that these extraneous factors can be discounted 
as discriminative because only ligand nature and concentration 
varied between stimulations. These built-in controls justify the 
simplifying assumptions from which the mathematical model was 
derived. 

The present findings certainly do not exclude the existence of a 
receptor-level hierarchy. Indeed, the co-existence of hierarchies at 
both the receptor and cellular levels may explain the data equally 
well. However, the hypothesis that the response hierarchy resides 
solely at the receptor level is contradicted by the present findings. A 
threshold hierarchy at the level of cellular activation would imply 
that the functional responses are located downstream from signal- 
ing pathways that share a common starting point (the triggered 
TCR or signalosome). Further detailed molecular studies of the 
CD3-complex phosphorylation cascade would be required to rule 
out a receptor-level hierarchy definitively. 

Our mathematical model elucidates the contrasting roles played 
by the dissociation rate and the affinity constant as determinants 
of T-cell functional sensitivity, with TCR and pMHCI molecular 
densities governing the limitation regime under which the kinetics 
operate. The disappearance of the optimum behavior as the sys- 
tem moves from the receptor-limited to the affinity-limited regime 
can be studied experimentally by increasing TCR density. This 
shift was first predicted by Van den Berg et al. (23) and confirmed 
experimentally by Gonzales et al. (41). Furthermore, the pattern 
of dependence on fc on and k 0 g shown in Figure 2B may account 
for the discrepancies between the findings of Kalergis et al. (42), 
who reported such an optimal dwell-time (= k~2), and the find- 
ings of Holler et al. (43), who reported monotone dependence. 
The shift along the fc on axis may also explain the observations of 
Irving et al. (26), who observed that affinity modulates the depen- 
dence on k 0 g. In the experimental system studied here, the TCR 
is stimulated under MHC-limited conditions. This means that 
the optimum-type dependence on k 0 g prevails. This is shown in 
Figure 6, where the calculated TCR triggering rate (at 10~ 5 M) is 
plotted against /c Q ff. The effect of k on is so slight that the curves 
overlap almost perfectly for the five ligands. Even under MHC- 
limited conditions, however, k on remains a crucial determinant 
of functional sensitivity, in view of the fact that the ratio of fc G ff 
and fc on determines the transition between affinity-limited and 



ligand-limited signaling, as illustrated schematically in Figure 7. 
Thus, the association rate and the dissociation rate play distinct 
roles in T-cell activation, consistent with reports of dominant 
effects associated with each of these biophysical parameters under 
different conditions (24, 44); these differences also explain earlier 
observations on the interplay of kinetic parameters and receptor 
numbers (41, 45, 46). 

In the presence of the wild-type MHCI/CD8 interaction, the 
ligands 3G and 3G8T are virtually optimal, whereas abrogation of 
MHCI/CD8 binding leads to a reduction in functional sensitivity. 
In virtue of the non-linear character of the curve, the magnitude 
of the reduction varies considerably across ligands. Two parame- 
ters (y on and j/r) suffice to capture these effects. The ILA ligand 
evokes a triggering rate approximately one order of magnitude 
less than optimal under C1R-A2null conditions. Nonetheless, the 
data show that strongly attenuated but discernible responses can 
be elicited at sufficiently high presentation levels, in keeping with 
Wooldridge et al. ( 14), who reported that a ligand could be physio- 
logically significant even at an estimated ~2.1 orders of magnitude 
below the optimal ligand for a responding T-cell clone. 

Although all ligands studied here become less potent when 
presented by C1R-A2null cells, the model predicts that this is 
not necessarily the case. Indeed, the functional sensitivity to a 
ligand marked "X" in the inset of Figure 6 is unaffected when 
the MHCI/CD8 interaction is abrogated. Moreover, ligand "Y" 
becomes optimal in the absence of the CD8 effect (so that the T-cell 
has a higher functional sensitivity to this ligand under CD8-null 
conditions) . The existence of such ligands would imply that T-cells 
can "tune" to distinct cognate ligands by up- or down-regulating 
their CD8 coreceptors. This effect would greatly amplify the ability 
of T-cells to reconcile extensive cross-reactivity with the avoidance 
of self-recognition. 

The hypothetical ligands "X" and "Y" are strongly hetero- 
clitic. Peptides of this nature occur markedly less frequently than 
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FIGURE 7 | Interplay between kinetic parameters and receptor 
densities. (A) Affinity-limited: theTCR triggering rate is proportional to the 
densities of both receptors. (B)TCR-limited: theTCR triggering rate is 
proportional to theTCR density and is monotone decreasing in k on . (C) 
MHC-limited: theTCR triggering rate is proportional to the MHC density 
and has an optimum k„„.Jhe affinity K D = k 0 „/k m governs the transitions 
between these kinetic regimes. 



homoclitic ligands, which makes them more challenging to detect 
and characterize. Nonetheless, work is ongoing to demonstrate the 
existence of such anomalous ligands, using combinatorial peptide 
library screening combined with importance sampling ( 14) in the 
presence of different extracellular contributions from the CD8 
coreceptor (29). 

The effect of CD8 on the value of the TCR triggering thresh- 
old is consistent with the conclusions of Van den Berg et al. (39) 
and may be due to the association of TCR/CD3 with protein tyro- 
sine kinases such as p56lck, which expedites the immunoglobulin 
family tyrosine-based activation motif phosphorylation sequence 
(47). This agrees with the observation that the CD8 aft hetero- 
dimer is more potent as a coreceptor than the aa homo-dimer 
(48, 49), perhaps due to interactions with other signalosome 
components mediated by the palmitoylated CD8 ft chain (50), 
which interacts with myristoylated p56lck. Docking of CD8 to 
the TCR/CD3 complex may also be driven by initial activation 
(51). The kinetics of CD8 binding after TCR/pMHC contact 
was analyzed previously by Van den Berg and Sewell (52). More 
recently, Mukhopadhyay et al. (53) developed a model that explic- 
itly accounts for the interplay between the CD3 phosphorylation 
cascade and p56lck, ZAP-70, and CD45. 

The dissociation rates determined by SPR were assumed to be 
indicative of the normal interaction in this study. In reality, these 
values might be more true to the case where the MHCI/CD8 inter- 
action is absent. Assuming the latter would not dramatically alter 
the qualitative conclusions of this study, but would lead to a slight 



adjustment in the interpretation of the parameter estimates. Most 
importantly, the value found for the TCR triggering threshold Tr 
would be adjusted down from 17.46 to 15.26 s. 

Values based on SPR experiments are best regarded as 3- 
dimensional, representative of ligands in solution. However, the 
TCR/pMHCI interaction takes place in a 2-dimensional envi- 
ronment, which essentially reduces spatial degrees of freedom of 
molecular motion and introduces dynamics related to the forces 
that constrain the molecules to this environment (54). Conse- 
quently, rate constants may be different in two dimensions as 
opposed to three dimensions. In particular, 2-dimensional dis- 
sociation rates can be substantially faster (55). The ratio between 
the 3-dimensional and 2-dimensional affinities is a length mea- 
sure, denoted h and called the confinement length (56, 57). Wu 
et al. demonstrated that h is proportional to the range of motion 
available to the free forms of the interacting ligands along the 
spatial axis perpendicular to the two parallel membranes. Thus, 
the intermembrane separation provides an upper bound, and if 
the ranges of motion do not differ too much, it is reasonable to 
assume that the confinement length is roughly the same for all 
mutants involved. 

The use of a fixed receptor duration might appear simplistic 
given the complexity of the events that are required to trigger a 
TCR/CD3 complex, but the law of large numbers provides indi- 
rect support for a simple assumption: the true waiting-time-till- 
triggering is a composite of a large number of stochastic durations 
which tends to regress it on the mean value. In fact, the number of 
terms is itself a random variate due to the possibility of alternative 
routes in the activation pathway of the complex (such processes 
are termed "compound stochastic") . The complexity illustrated in 
Figure 1 effectively underpins a simple assumption, viz. that there 
is a fixed, stereotypical receptor threshold Tr. 

The parameters Tr and w act are both subject to further modula- 
tion. By varying the levels of the CD8 coreceptor, surface molecules 
such as CD45, and cytosolic concentrations of kinases, phospho- 
rylases, and linker proteins in the immediate vicinity of the CD3 
complex, the T-cell can adjust the effective value of the recep- 
tor triggering threshold T R (29, 30, 38, 39, 58-61). The value of 
the cellular activation threshold w act is modulated by costimu- 
lation ("signal 2"). Moreover, the T-cell can actively modulate 
its receptiveness to this signal by adjusting the relative levels of 
CD28 and CD152 (52, 62). Furthermore, the maintenance of tol- 
erance requires continual dynamic tuning of the cellular activation 
threshold (63, 64). Conceivably, these mechanisms provide an 
additional "multiplier" that modulates w act over and above the 
factors that set w act to different values for different responses. 

The data on intercellular response variability (Figure 3) sup- 
port the assumption that the distribution of responsiveness in 
terms of MFI can be attributed to intercellular variation in the 
triggering threshold. However, the observed variability may also 
arise from differences in the numbers of TCR and pMHCI mole- 
cules present in the T-cell/ APC interaction area. Mathematically, 
this model is slightly more involved because interactions can now 
be distributed over different kinetic regimes. Nonetheless, apart 
from this technical difficulty, the resulting model is mathemati- 
cally equivalent; this follows from the properties of the log-normal 
distribution. Thus, the analysis is not substantially affected by this 
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alternative explanation of the intercellular variability. Experimen- 
tally, the two hypotheses could be distinguished by following an 
individual T-cell over a series of interactions and responses. If cel- 
lular thresholds for the various responses remain largely constant, 
there would be little event-to-event variation according to the for- 
mer hypothesis; in contrast, the receptor copy number hypothesis 
would predict considerable event- to -event variation. 

In addition, there may be variability with respect to antigen 
presentation levels across the experimental APC population, due 
to differential antigen exposure and stochastic effects in anti- 
gen acceptance from the incubation medium. However, it can 
be shown that the resulting variance is not such as to make 
a substantial contribution to overall variability, for instance by 
following the argument presented in Appendix C of Van den 
Berg et al. (22). 

When the TCR/pMHCI ternary complex dissociates, the cova- 
lent modifications and multimeric aggregations that the CD3 
complex undergoes during the gradual transition to signalosome 
status are reversed, both by thermal agitation and the action of 
phosphorylases. In the present model, it has been assumed that this 
"reset" occurs much more rapidly than the typical time required 
for the TCR/CD3 complex to encounter the next non-null pMHCI 
ligand. However, given that the molecules have to drift away from 
each other through diffusion in the 2-dimensional arena of the T- 
cell/APC interface, there is a possibility of rapid rebinding to the 
same pMHCI molecule before the CD3 complex has had sufficient 
time to "reset." This rapid rebinding is equivalent to a diminish- 
ment of the effective off-rate; a basic model capturing this effect 
has been proposed by Aleksic et al. (65). Extending the present 
model with this effect is straightforward in principle, albeit at the 
cost of an additional parameter. Furthermore, the number of com- 
plexes with other pMHCI species is ignored in equation (1). This 
is valid if the vast majority of all other ligands consists of null 
agonists. An analysis that takes the entire presentation profile (all 
pMHCI species present) into account can be found in Van den 
Berg etal. (22,23). 

A further aspect of the cellular activation threshold that has 
been left implicit in the present treatment is the duration of the 
T-cell/ APC interaction. If this contact is initiated at to and is ter- 
minated at fi, the cumulative signal transmitted to the cytosolic 



signaling machinery is given by an integral: 

Q = I ' W(x)dx. (16) 

J t 0 

In fact, this quantity Q is presumed to be directly compared to 
a threshold. However, Q/( 1 1 — to) is proportional to the activation 
threshold as defined in the present study. Indeed, variability in the 
duration of the contact is one of the sources of stochasticity under- 
lying the log-normal response curve used to fit the dose-response 
curves. 

A related phenomenon is that of TCR down-regulation (66- 
68), whereby triggered TCRs are gradually removed from the cell 
surface. This affects the integral in equation (16); Van den Berg et 
al. (23) discuss how the model can be extended to take this effect 
into account. TCR down-regulation can in fact be exploited by the 
T-cell to gauge both fc on and k a ff independently. As indicated in 
Figure 2A, if the system starts in the MHC-limited regime, grad- 
ual removal of TCR molecules will result in a transition to the 
TCR-limited regime, where a sudden drop of the instantaneous 
signal results. At the point of change-over, Rj^ Zj. This means 
that the T-cell can, in principle, glean a rough estimate of the bio- 
physical rate parameters k 0 ff and k on by combining information 
from the gross signal and the number of down-regulated TCRs 
at the transition point. This mechanism can only function in the 
ligand-limited regime, not in the affinity-limited regime. 

In conclusion, the present study confirms that a response hier- 
archy exists with respect to the strength of TCR stimulation 
required to elicit various cellular responses. Moreover, a combi- 
nation of experimentation and mathematical modeling indicates 
that this hierarchy resides at the cellular level rather than at the 
level of the individual receptor molecules. Functional sensitivity 
is generally enhanced by the CD8 coreceptor, although the the- 
ory indicates that the MHCI/CD8 interaction can depress TCR 
signaling for certain heteroclitic ligands. 
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FIGURE Al | Kinetic measurements of ILA1 TCR interactions with the 
agonist pMHCI ligand 3G8T. SPR kinetic measurements were performed 
with five serial twofold dilutions of the ILA1 TCR, starting at 50 /iM, flowed 
over immobilized 3G8T/HLA-A2 complexes. Injections were performed in 
duplicate for each TCR concentration. The plots show data (thin lines) and 
curve fitting (thick lines) representative of two separate experiments. 
First-order exponential association rate (k m ) and dissociation rate (Ar off ) were 
fitted simultaneously by least squares. Values for k 0 „, k 0 „, and K D are shown 
in Table 1. 
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